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abstract 

In this paper we develop an approach to conformal geometry of piecewise flat metrics 
on manifolds. In particular, we formulate the combinatorial Yamabe problem for piecewise 
flat metrics. In the case of surfaces, we define the combinatorial Yamabe flow on the 
space of all piecewise flat metrics associated to a triangulated surface. We show that the 
flow either develops removable singularities or converges exponentially fast to a constant 
combinatorial curvature metric. If the singularity develops, we show that the singularity 
is always removable by a surgery procedure on the triangulation. We conjecture that after 
finitely many such surgery changes on the triangulation, the flow converges to the constant 
combinatorial curvature metric as time approaches infinity. 

§1. Introduction 

1.1. In this paper we consider the class of piecewise flat metrics on a triangulated 
manifold and develop an approach to piecewise linear (PL) conformal geometry of these 
metrics. We also formulate the combinatorial Yamabe problem for piecewise flat metrics. 
The notion of conformal class of a metric was originated from Riemannian metrics. Given a 
Riemannian metric gij on a smooth closed manifold, the conformal change of the metric is 
a new Riemannian metric of the form ugij where u is some smooth positive function defined 
on the manifold. The Yamabe problem states that there exists a constant scalar curvature 
metric in the conformal class of any Riemannian metric. This was solved affirmatively 
for closed Riemannian manifolds in [Ya], [Tu], [Au] and [Sc]. The solution was considered 
as a milestone in application of non-linear partial differential equations in geometry. In 
trying to develop the analogous PL conformal geometry, one faces the task of defining 
the notion of PL conformal class of a PL metric. The most naive approach is to take 
a local (infinitesimal) point of view. Namely, since the conformality in smooth category 
means infinitesimal invariance of angles, thus PL conformal geometry should preserve the 
of measurment of angles in the PL metrics even at the singularities. However, we are not 
able to come up with any reasonable results from this point of view. It is probably due 
to the that fact that local geometry of smooth metrics is governed by linear algebra and 
the local geometry of PL metric is governed by the geometry of simplexes. We propose to 
approach the analogous PL conformal geometry of PL metrics from global point of view 
in this paper. Namely, the conformal class of a smooth metric comes from the action of 
smooth functions on smooth metrics. Thus the PL conformal class of a PL metric should 
be a class of PL metrics related by the action of functions defined on the set of all vertices. 
Our main results indicate that the PL conformal geometry captures some of the main 
features of the smooth conformal geometry on surfaces. The situation in PL theory is 
more complicated in the sense that there are triangulations on surfaces so that any PL 
metric associated to the triangulatinos has none constant PL scalar curvature. This is 
mainly caused by the fact that a PL metric is supported by many different triangulations. 
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In the case of surfaces, we prove that the existence of a constant PL scalar curvature metric 
is a purely topological condition on the triangulation (theorem 1.1). We propose to find 
the constant PL scalar curvature metric in the PL conformal class by a system of ordinary 
differential equations. These equations seems to be the right analogy of the Yamabe flow. 
We establish some of the basic properties of the flow. Especially, we prove that the flow 
will not develop essential singularities in finite time and converges exponentially fast to the 
constant PL curvature metric if no singularity develops. The most interesting property of 
the flow is that when the flow develops singularity, the flow naturally suggests a way to do 
surgery on the underlying triangulation. This indicates that the flow tends to search for 
both the best piecewise flat metric and the underlying triangulation. 

1.2. We begin with a dictionary between piecewise linear (PL) and smooth theory. In 
the smooth theory, we start with a smooth manifold N n and whereas in the PL theory, 
we start with a triangulated manifold (M n ,T) where T is the triangulation. The natural 
analogy of the smooth functions on N n are the functions defined on the set of all vertices 
V in the triangulation T. Similarly, the analogy of tensors on N n are functions defined 
on the set of all z-simplexes in the triangulation. For instance, the de Rham cohomology 
and simplicial cohomology are such examples. We define a PL metric associated to the 
triangulation to be a positive real valued function defined on the set E of all edges (i.e., 
1-simplexes) of T so that for each n-simplex in T, the restriction of the function to the 
1-simplexes of the n-simplex is the edge lengths of some Euclidean n-simplex. The last 
condition is the realizability condition which is given by the Cayley-Menger matrices for 
Euclidean n-simplexes. We may think of a PL metric as a metric on the manifold so that 
its restriction to each simplex is isometric to a Euclidean simplex. Just like the space of 
all Riemannian metrics on a smooth manifold is convex, it can be shown that the space of 
all PL metrics associated to a triangulation is a convex set. Indeed, it is an observation 
of Igor Rivin that the set {d 2 : E — > R| where d : E — > R>o is a PL metric associated 
to T} is convex in Ti E . For any Riemannian metric, the curvature tensor assigns to every 
four tangent vectors at a point a number. For the PL metric, there are many notions 
of curvatures (see [CMS], [CE] and others). The most natural curvature seems to be the 
function K which assigns to each (n-2)-simplex a n _2 the number K(a n -2) = 2n — a where 
a is the sum of all dihedral angles of n-simplexes at the (n-2)-simplex. It is shown in [St] 
that if the K(a n -2) > for all (n-2)-simplexes, then the manifold has finite fundamental 
group. This is the analogy of the Myers' theorem in the smooth theory. The combinatorial 
scalar curvature (or PL scalar curvature) of a PL metric d is a function Sd : V — > R 
whose value at the vertex v is Sd(v) = J2 tTn _ 2>v K{a n -2)vol{a n -2) where the summation 
is over all (n — 2)-simplexes having v as a vertex. If n = 2, then, as a convention, we set 
vol(v) = 1. Note that if A is a positive number then S\d = \ n ~ 2 Sd and also that the total 
scalar curvature J2 veT o S(v) is exactly Regge's proposed approximation to the Einstein 
action (see [Re]). Furthermore, the PL scalar curvature Sd depends not only on the metric 
d (as a metric in point set topology) but also on the choice of the underlying triangulation 
T. If the manifold is a surface, then the scalar curvature S(v) is the curvature K(v). 

Using this dictionary, it is now natural to define a PL conformal factor for a PL 
metric to be a positive function u defined on the set V of all vertices. Given a PL metric 
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d : E — > R>o, we define the PL conformal change of d to be the new metric u*d : E — > R >0 
where -u * d(vv') = u(v)u(v')d(vv') and uu' is the edge with vertices i> and ?/ so that u* d 
is realizable on each n-simplex. The set of all PL metrics of the form u * d for a fixed d is 
called the PL conformal class of the PL metric d. 

The combinatorial Yamabe problem asks if there is a constant PL scalar curvature 
metric within the PL conformal class of any PL metric. As we will see, there are topological 
obstructions for the existence of the constant scalar curvature PL metrics associated to a 
triangulation even in dimension 2 (see theorem 1.1). The goal is to establish the existence 
of constant PL scalar curvature metric when the topological obstructions vanish. 

1.3. In the case of a triangulated surface, we propose to approach the combinatorial 
Yamabe problem using a system of ordinary differential equations. Here is the setup. Let 
T be a triangulation of a closed surface M and V and E be the set of all vertices and 
edges in T. A PL metric associated to T is a positive function d : E — > R so that it 
satisfies the triangular inequalities on the three edges of any triangle in T. Fix the PL 
metric d on the surface. Let u : V — > R>o be a conformal factor. For simplicity, let us 
write V = {vi, ...,vn} and Ui = u(vi). If ViVj is an edge in T, then d(viVj) is denoted 
by dij. We define the combinatorial Yamabe flow to be the following system of ordinary 
differential equations. 



(1.1) = 

ut(0) = 1 

Here the curvature Ki is the K curvature of the metric u*d at time t at the i-th vertex Vi. 

The Gauss-Bonnet theorem for PL metrics states that J2iLi = ^xiM). Thus the 
average curvature K av is 2nx(M)/N . In particular, the constant combinatorial curvature 
PL metric has curvature K av at all vertices. The combinatorial Yamabe problem asks if 
this metric exists in the PL conformal class. 

There is a combinatorial obstruction for the existence of the constant curvature PL 
metric associated to a triangulation T. In fact the following holds. (For a finite set X, we 
use \X\ to denote the number of elements in X.) 

Theorem 1.1. Fix a triangulation T of a closed topological surface M . There exists 
a constant PL scalar curvature metric associated to T if and only if for any proper subset 
I of the vertices V of T , 

(*) \Fi\/\I\ > \F\/\V\ 

where F is the set of all triangles in T and Fj is the set of all triangles having a vertex in 
I. 
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Furthermore, the condition (*) holds for all triangulations of surfaces of non-negative 
Euler characteristic. 

For simplicity, we call a triangulation which supports a constant curvature PL metric 
admissible. 

There are none admissible triangulations on each closed surface of negative Euler char- 
acteristic. It can be shown that given any triangulation of a surface, there is a subdivision 
of it which is admissible. However, we are not sure if the subvision can be choosen to be 
an iterated barycentric subdivision or even the 2-dimensional subdivision which replace 
each triangle by four triangles formed by the barycenters of the vertices and edges. On the 
other hand, if our conjecture in section 1.3 holds, then every triangulation of the surface 
is equivalent to an admissible one having the same number of vertices. 

Our main results concerning the combinatorial Yamabe flow can be summarized in 
the following theorems. 

Theorem 1.2. Fix a piecewise flat metric d on a triangulated surface (M,T). The 
following holds. 

(a). Under the combinatorial Yamabe flow, the curvature evolves according to a com- 
binatorial heat equation of the form 



(i-2) ^r = £ c ^ 

where the matrix [cy]jvxiv is symmetric and semi-negative definite. In particular, the total 
curvature Ki{t) 2 is decreasing in time t. 

(b) . The combinatorial Yamabe flow is variational. If we change the variable U{ to 
Wi = log then the combinatorial Yamabe flow is the negative gradient flow of a locally 
convex function in w = (wi, wn)- 

(c) . (local rigidity) The curvature map K : {u = (u\, ...,un) G R>ol Yli=i u % = 1? u *d 
is again a PL metric } — > {(ki, fejv)| Z)£Li = 27i~x(M)} sending a conformal factor u 
to the curvature of the metric u * d is a local homeomorphism. 

We do not know if this local rigidity can be improved to be a global rigidity result. 

To understand the long time behavior of the flow, we introduce the normalized equa- 
tion 



(1.3) 



~dt 



-Ui(Ki - K av ) 



«i(0) = 1 



4 



The two equations (1.1) and (1.3) are equivalent in the sense that Ui(t) solves (1.1) 
if and only if e 27rx ^ M ) t / Ai Ui(t) solves (1.3). Furthermore, if u{t) = (iti(t), itjv(t)) is a 
solution to the normalized equation (1.3) then the product Y[iLi u i(t) = 1 f° r an time t. 
To understand the asymptotic behavior of the solution of (1.3), we have to analysis the 
potential formation of singularities in the normalized equation (1.3). There are only two 
types of singularities which may occur. Suppose u(t) is a solution of (1.3) in the time 
interval [0, L) where L < oo. We say the solution develops an essential singularity at time 
L if there is an index i and a sequence of time t n approaching L so that limt n ^L Ui(t n ) = 0. 
We say the solution develops removable singularity at time L if there is a sequence of time 
t n approaching L so that Ui(t n ) remains in a compact set in R >0 for all % and there is a 
triangle AviVjVk in T which degenerates into a line segment as t n — > L. The last condition 
means that the triangular inequality UidikUk + UkdkjUj > UidijUj becomes equality as 
t n — > L for some indices k. 

Theorem 1.3. (a) For any triangulated surface, the normalized combinatorial Yam- 
abe flow will not develop essential singularity in finite time. 

(a) If the triangulation is admissible, the normalized combinatorial Yamabe flow will 
not develop essential singularity at time infinity. 

If removable singularity occurs, say the vertex Vk is moving toward the interior of the 
edge ViVj, then we change the original triangulation T as follows. Let v\ be the fourth 
vertex in T so that ViVjVi forms a triangle in T. Then we replace the triangulation T by 
a new triangulation T' which is obtained from T by deleting the edge v^j and adding 
a new edge v^vi. In particular, two triangles AviVjVk and AviVjVi in T are replaced by 
two new triangles Av^viVi and AvkViVj. This is the most common surgery operation on 
triangulations of surfaces. After the combinatorial surgery on T, we run the combinatorial 
Yamabe flow on T" with initial metric coming from the final stage of metric on T at time 
L. We conjecture that after finitely many such surgery operations on the triangulation 
T, the combinatorial Yamabe flow converges exponentially fast to the constant curvature 
metric. This is supported by the following. 

Theorem 1.4. If no singularity develops in the normalized combinatorial Yamabe 
flow, then the solution converges exponentially fast to a constant curvature PL metric as 
time approaches infinity. 

Under the normalized flow, the distance function restricted to the 1-skeleton of the 
triangulation will stay in the same quasi-isometric class in any compact time interval in 
R>o- The geometry of the quasi-isometry class changes only when the essentialy singulariy 
develops at time infinity. 

If we start with a single Euclidean triangle, and define the curvature Ki at the i-th 
vertex to be tx — 9i where 6i is the inner angle at the i-th vertex, then it can be shown 
that the corresponding normalized equation (1.3) has solution for all time and converges 
exponentially fast to the equilateral triangle as time tends to infinity. 

The following is likely to hold for the combinatorial Yamabe flows on surfaces. 
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Conjecture, (a) A constant PL scalar curvature metric is unique in its PL conformal 
class. 

(b) The combinatorial Yamabe flow will not develop essential singularities on any 
triangulation. 

(c) The combinatorial Yamabe flow will converge to the constant curvature PL metric 
after finite number of surgeries on the triangulation. 

There are some evidences indicating that the combinatorial Yamabe flow will develop 
removable singularities for some PL metrics on admissible triangulations. But we do not 
have the explicit examples yet. 

It seems there is an interesting similarity between the singularity formation in combi- 
natorial Yamabe flow on surfaces and Hamilton's Ricci flow program in dimension 3. 

1.4. There are several interesting questions concerning the combinatorial Yamabe 
problem in higher dimension. To carry out the same program in higher dimension, the key 
ingredient that is missing is the following local rigidity property. 

Question. Given an Euclidean n-simplex a n and an (n-2)-face o~ n -2 of a n , let 
o(cr n _ 2 ,a n ) be the dihedral angle of the n-simplex at the (n-2)- simplex. Now suppose 
d : E — > R>o is the edge lengths of the n-simplex. Consider conformal factors u : V — > R>o 
and the new PL metric u * d on the n-simplex. Define the function S u *d '■ V R n to be 
S u *d(v) = J2 a - ri 2>v a(a n _2, o~ n )v ol(a n -2) in the metric u*d. Then for fixed metric d, S u *d 
is a smooth function of u. Is the rank of the Jacobian matrix of the map S u *d considered 
as a function defined in an open set in H v to H v always equal to n — 1 ? 

An affirmative answer to this question for n > 3 will give a strong evidence that the 
higher dimensional combinatorial Yamabe flow exists, i.e., one would call equations of the 
type = f( u i)S( v i) the combinatorial Yamabe flow where f(x) is some universal no- 
where zero function depending on the dimension n. The function / is so chosen that the 
corresponding matrix [^-f(uj)] is symmetric. This question is also related to the follow- 
ing. Given a PL metric d associated to a triangulation, consider the set of all conformal 
factors u for d. Let R(u*d) = J2 a _ 2 K{a n -2)vol{a n -2) be the Regge action of the metric 
u* d. Is the Hessian of R(u * d) considered as a function of u semi-positive definite? 

The other related question we are considering now is the combinatorial Yamabe flow 
in hyperbolic or spherical background metric, i.e., we use hyperbolic or spherical simplexes 
instead of the Euclidean ones. The work of [CL] suggests that in general it is easier to 
work in hyperbolic geometry. 

There are many interesting questions in the field of PL metric theory. For instance 
what is the right analog of the Ricci curvature? What is the Laplacian operator on the space 
of all functions defined on the i-simplexes? Is it true that positive PL scalar curvature 
carries topological information? To be more precise, suppose M is a 3-manifold which 
supports a positive PL scalar curvature. Is it ture that the fundamental group of M 
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contains no none-abelian surface group? In the smooth case, the result was proved by 
Schoen and Yau [SY] . 

1.5. In [CR], Cooper and Rivin define the concept of combinatorial scalar curvature 
differently. In their definition, the combinatorial scalar curvature is invariant under the 
scaling of metric. On the other hand, the scalar curvature in Riemannian geometry has 
the property that S(kgij) = l/kS(gij) for any positive constant k. (In our case, the PL 
scalar curvature satisfies S(kg) = k n ~ 2 S(g).) It can be shown every topological 3-manifold 
has a triangulation and an associated PL metric whose PL scalar curvature is positive in 
the sense of [CR] . Thus the positivity of PL scalar curvature defined in the sense of [CR] 
does not carry any topological information. However, at each vertex, the scalar curvature 
defined in [CR] is a measurment of the infinitesimal rate of change of the volume of a small 
ball centered at the vertex. Thus, it is a combinatorial analogy of the scalar curvature 
in the infinittesimal sense. The work [Ga] follows the approach of [CR] and defines a 
combinatorial Yamabe flow for ball packing metrics. The PL scalar curvature defined in 
this paper is more closely related to Regge's calculus. 

1.6. The paper is organized as follows. In §2, we establish some basic properties of the 
geometry of triangles. In §3, §4, and S5, we prove 1.1, 1.2, 1.3, and 1.4. In the appendix 
A, we provide a detail calculation involved in the proof of theorem 2.1. In appendix B, we 
use the feasible flow theorem for network flow to establish the necessity part of theorem 
1.1. 

1.6. This work is supported in part by the NSF. We thank Ben Chow, X-S. Lin and 
D. Sullivan for discussions. We thank the referee for making some nice suggestions. 

§2. Geometry of Euclidean Triangles 

For simplicity, let A = {(xi, £2, £3) G R^ol 2 -* + x j > x k, where i, j, k are distinct} 
be the set of points whose coordinates satisfy the triangular inequalities. We use R>o to 
denote the set of all positive numbers, and Av\V2V^ to denote the triangle having vertices 
vi, V2, and V3. If X is a finite set, we use \X\ to denote the number of elements in X. 

2.1. Fix three numbers di, dj, dk satisfying triangular inequalities, i.e., (di, dj, du) G A 
and three vertices Vi,vj,Vk- Consider the conformal factor u : {vi,vj,Vk} — > R>o so that 
(di/ui,dj/uj,dk/v,k) G A where u r = u(v r ). Define x r = d r u s ut where {r, s, t} = {i,j,k}. 
Due to the choice of w r 's, (xi, xj, Xk) G A. Construct a Euclidean triangle AviVjVk so that 
the length of the edge opposite to v r is x r . Let 9 r be the inner angle at the vertex v r . 
Then 9 r is a smooth function of (ui, uj, u^)- 

Theorem 2.1. The 3x3 matrix [^ !L u s ]sx3 is symmetric, semi-negative definite and 
has rank 2 whose null space is {(t,t,t) G R 3 | t G R}. Furthermore, if we let z. and 

a rs = —x r x s cos(6> t ) where {r, s, t] = k} and let A be the area of the triangle AviVjVk, 
then [|£ w s]3x3 = -2^[a rs ] 3x3 . 
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The proof is a straight forward computation. We defer it to the appendix A. It is easy 
to see that the matrix [a rs ]3 X 3 is symmetric and semi-positive definite. Indeed, since Xi = 
Xj cos 9k+Xk cos Oj, we see that the sum of entries in every row is zero. Thus <iet([a rs ]3 X 3) = 
0. On the other hand, the diagonal entries are positive so are the determinants of the 
principal 2x2 submatrices. Thus the 3x3 matrix [a rs ]3 X 3 is semi-positive definite whose 
null space is the diagonal in R 3 . 

Corollary 2.2. The matrix [^f L u s ]3x3 depends only on the three inner angles 9i,9j 
and 9k of the triangle. In particular, if the conformal factor u = (ui,Uj,Uk) varies in a 
region so that the inner angles of the triangle AviVjVk lie in a compact set in the open 
interval (0, tt), then there is a positive constant X so that the negative eigenvalues of the 
matrix [^ 1 u s ]s x s are less than —A for all u in the region. 

2.2. Using the same notations as above, let w r = logu r for r = k. Then the 
3x3 matrix [§^-]3 X 3 is symmetric and semi- negative definite of rank 2 whose null space 
is {(£, t, t) G R 3 |t G R}. Consider the non-convex space W = {w = (wi,Wj,Wk) G 
R \(die~ Wi , dje~ Wi , dke~ Wk ) G A}. The space W is simply connected since it is the image 
of the convex space {(n, r,-, r^) G R 3 | (diVi, djrj, dkVk) G A} under the homeomorphism 
h(ri,rj,rk) = (— logr^, — logrj, — logr/j). Furthermore, if w G W, then w + (t, t, t) is again 
in W for any real number t. Take any vector (a$, aj, ak) G R and let O = (a^ — 9i)du)i + 
(cij — 9j)dvjj + (afc — 9k)dwk be a smooth 1-form defined on the space W. It is closed due 
to the symmetry of the matrix [§^-]3x3- Define a smooth function F(w) : W — > R by 



This is well defined since the form O is closed and the space W is simply connected. 

Corollary 2.3. The function F : W — > R is a locally convex function so that it 
becomes a locally strictly convex function when restricted to the planes {(wi,Wj,Wk) G 
+ Wj + Wk is a constant }. Furthermore, 

(a) , if di + aj + ak = 7T, then F(w + (t,t,t)) = F(w) for any real number t. 

(b) . if (fli, aj, ak) G R^q and ai + aj + ak = tt, then for any sequence of points in 



W, limsup n ^ 00 F(tt;( n )) = oo if and only if limsup n ^ oc max {rjS}c{i)J)fc} (|«4. nJ = 



Proof. The local convexity is due to the fact the Hessian of F is the matrix — [J^]3 X 3 
which is semi-positive definite. Furthermore, since the matrix is positive definite on the 
linear space Wi + Wj +Wk=0, it follows that F is locally strictly convexity on the planes. 

To see (a), it suffices to verify that Jj^+( d ' d ' d ) q _ q_ Now take the line segment 

tw + (1 — t)(d, d, d) in W to evaluate the line integral q -\y e g nc [ integral 

becomes f Q ((a^ + aj + ak) — {9i + 9j + 9k))dt = since the integrant is 0. 

To see (b), let IT : R 3 — > {(wi, W2, w^) G H 3 \wi + w 2 + W3 = 0} be the orthogonal 
projection. Then due to the condition ai + aj + ak = 7r and part (a), we have F(w) = 



(2.1) 




00. 
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F(H(w)). Also, we have limsup n ^ 00 max{ r s }. C {j ^/ c }(]tUr — Ws |) = oo if and only if their 

projections n(W( n )) is unbounded. Thus if limsup n _ 00 m.ax^ s y C [ ijjjk y (\wP-wi n) \) is 
finite, then the function F(w n ) remains bounded. 

Now suppose limsup n ^ QO rnax{ T . jS } C { ijJ - j fc} — wi™' | = oo. By taking a subsequence 
if necessary, we may assume without loss of generality that lim n ^ 00 (u^ n ' ) — Wj) = oo. By 
part (a), we may further assume after adding by a diagonal vector (d ni d n ,d n ) that 
w (n) j s m ^ ne hypersurface {{wi,Wj,Wk) £ VF|e _u;i + e _Wj + e _Wfe = 1}. In particular, 

(n) (n) — — — 

Wr > 0. Thus, lim n ^ 00 = oo. By the constraints e Wi + e + e Wfc =1 and 
the triangular inequalities die~ Wi + d r e~ Wr > d s e~ Wa for {r, s} = {j, k}, both and 
wj^ are bounded from above. This shows that the triangle AviVjVk with edge lengths 

(n.) - (tl) (ri) (n) (n) , (n) 

c^e™-? fe , (ij-e^i +u; fc and d] i e Wi i degenerates into a half-line since two of the 
edge lengths tend to infinity and the third remains bounded. Thus 9i{w^) tends to 

w (n) li/™" 1 

zero. Therefore the dominate term in the integration J Q O is J Q i (a^ — 9i)dwi since 
the other two integrals are bounded. Because ai > and lim n 6i(w^) = 0, we see that 
lim^oo F(w^) = oo. QED 

2.3. Remarks. The situation in PL conformal geometry is very similar to the ap- 
proach to circle packing by Colin de Verdiere in [Cv] . The complication in the combinatorial 
Yamabe problem is caused by the fact that the space W is not convex and the function 
F : W — > R is not proper. 

§3. A Proof of Theorem 1.1 

That the condition (*) is sufficient was established in [CL] theorem 1.1 where we take 
the weight function = 0. Indeed, condition (1.3) for <fi = in theorem 1.1 in [CL] 
is exactly the condition (*). Thus in this case, there is a very special type PL metric 
associated to the triangulation whose curvature at each vertex is 2nx(M)/\V\. The metric 
is obtained by assigning to each vertex vi a positive number n and define the length at 
the edge ViVj to be + r j . 

To show that the condition is also necessary, we use the feasible flow theorem for 
network flow. Since the method of the proof is irrelevant to the rest of the paper, we defer 
the proof to appendix B. 

The fact that the condition (*) holds for all triangulations of surfaces of non-negative 
Euler characteristic was established in [CL]. 

§4. Proofs of Theorems 1.2. and 1.3 

4.1. To derive the evolution of the curvature Ki for the combinatorial Yamabe flow 
= —UiKi, we note that the evolution of an individual inner angle, say 9{ k in a triangle 
AviVjVk at the vertex Vi, is the following. 

dffj k __ d6j k du r 

dt ^ du r dt 

r=i,j,k 
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80 jk 

(4.i) =- y: ^~u r K r 



8u r 

r=i,j,k 



Now the curvature Ki at the i-th vertex is 2n — J2j k ®l where the summation is over 

dKj - ■">■" 

dt 

the curvature evolution equation is 

N 



all triangle AviVjVk in T. Thus = — J2j k ~ir- By (4-1) an d theorem 2.1, we see that 



df ~ 53 53 "' A '' " 53 r " A '' 
We can write the coefficient matrix as a sum of matrices, 

[Cir]NxN = 53^ 7^ u r]NxN 
3,k r 



jk 



where the sum is over all ordered edges {j, k}. The matrix A 3, = [-q^-u^n^n is symmet- 
ric and semi-positive definite with only one 3x3 submatrix being non-zero and all other 
entries being zero by theorem 2.1. Thus the matrix [c^] is symmetric and semi-positive 
definite. Furthermore, by the same identity above and theorem 2.1, we see that the rank 
of [cij] is N — 1 whose null space is spanned by the vector [1, 1, 1]. 
Take the time derivative of G(t) = Yl!i=i We get 

dC N 

— = 2j2K i K' i = 2j2c ij K i K j <0. 

i=l i,j=l 

Thus the total curvature square is decreasing along the solution. This establishes part (a) 
of theorem 1.2. 



4.2. To show part (b), if we change the variable Ui to Wi = logiti, then the combinato- 
rial Yamabe flow becomes = —Ki. Note that the symmetry of the matrix [^^Ug^xs 

implies that the N x N matrix [^ 1 }nxN is symmetric and is semi-positive definite. Its 
rank is (N-l) and its null space is {(t, t, t)\t E R}. 

Define W to be the space {w = (wi, ..., wn) £ H N \ if ViVjVk forms a triangle in T, 
then (wi, Wj, Wk) G W}. Note that W is not a convex space but is simply connected since 
it is the image of the convex set {(ri, r/v) e R>ol if ViVjVk forms a triangle in T, then 
diTi +djTj > dkTk } under the diffeomorphism (n, ...,rjv) — > (— logri, — logr/v). Also, 
if w e W, then w + (t, t, t) is still in W. 

Define a 1-form O = X^Li Kidwi on the space W. This form is closed since the matrix 
[^^jvxiv is symmetric. Thus the function 

f(w) = n 

Jo 
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is well defined on the simply connected space W. The function F is locally convex since 
its Hessian is the semi-positive definite matrix [|^]jvxAf and is locally strictly convex 
on the planes W C\ {w\wi + ... + wn = constant } since the matrix is positive definite 
when restricted to the plane + ... + wn = 0}. Due to ^ = Ki, it follows that the 

combinatorial Yamabe flow = —Ki is the negative gradient flow of a locally convex 
function defined on W. 

4.3. To see the local rigidity, we note that under the transformation from u to w, the 
space -u-space becomes the (N-l)-dimensional smooth manifold P = Wn{wi+...+WN = 0}. 
The curvature map K becomes the map sending w E P to the gradient of a locally strictly 
convex function F\ P defined on P. Due to the local strictly convexity, the gradient map is 
locally injective. Since the image of K is {(k\, fejv)| S^Li ^ = 2ttx(M)} which is also a 
manifold of dimension (N-l), by the invariance of domain theorem, the curvature map K 
must be a local homeomorphism. This establishes part (c). 

4.4. Remark. It can be shown that the curvature map K is a local diffeomorphism. 
However, we are not able to show that is globally a diffeomorphism. 

4.5. To prove part (a) of theorem 1.3, note that if Ui(t) solves the combinatorial 
Yamabe flow, then due to the universal curvature bound that (2 — \E\)n < Ki < 2n, we 
have 

l/ce~ ct < Ui (t) < ce ct 

for some positive constant c on the time interval [0, L) where the solution exists. Thus no 
essential singularity develops on any triangulation in finite time. 

4.6. To prove part (b) of theorem 1.3, since the triangulation is admissible, we find a 
PL flat metric associated to the triangulation so that its curvature is constant 2itx(M)/N. 
Let a{ k be the inner angle at vertex Vi in triangle AviVjVk in the metric. For the normalized 
equation u\ = —Ui(Ki — K av ), we make a change of variable Wi = logttj as before. By the 
normalization condition that Yli=i u i(t) = 1 for all time, we have Y^7=i w i(t) = 0- By the 
same calculation, we still have the fact that the N x N matrix [ d ^ Kr ^ w Kav ^ \NxN is again 
symmetric, semi-positive definite of rank (N-l) whose null space is 

R}. 

Now consider the space W as in subsection 4.2. We construct a specific 1-form O as 

N 

(4.2) o = £(]>>f - ei k )d Wi ) 

i=l j,k 

where 6j k is the inner angle of the metric u*d at the vertex Vi inside the triangle Av^jv^- 
By the choice of the inner angles aj k \ the 1-form O is exactly — K av )dwi. 

Let us now consider the associated function F(w) = f Q O defined on W. This func- 
tion is well defined since W is simply connected and O is closed. Furthermore, by the 
construction, the combinatorial Yamabe flow w'i = —(Ki — K av ) is the negative gradient 
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flow of F(w). In particular, if w(t) is a solution to the combinatorial Yamabe flow, the 
function h(t) = F(w(t)) is decreasing in time t. 

On the other and, we can rewrite the summation in (4.2) as follows, 

(af - e{ k )d Wi + (af - ef)dwj + (aj? - 0$)dw k 

L±ViVjV k kzl ~ 

where the sum is over all triangles in T. 

Since aj +a kl +a 1 ^ = re over each triangle, by corollary 2.3 (a), F(w+(d, d)) = F(w) 
for all w G W. Now suppose the normalized combinatorial Yamabe flow w\ = —(Ki — K av ) 
develops an essential singularity at time t = L < oo. Then, due to the normalization 
equation J^Li Wi(t) = 0, there is a sequence of time t n — > L and two indices i,j so that 
lim^^oo Wi(t n ) = oo and lim tn _ >00 Wj(t n ) = -oo. 

We claim in this case \im tn ^, L F(w(t n )) = oo. This will contradict the fact that 
F(w(t)) is decreasing in time t. 

To see that claim, we use the same argument used in the proof of corollary 2.3(b). 
Indeed, after adding w(t n ) by a vector of the form (d n , ...d n ), we may assume that w(t n ) 
is in the subspace {(wi, ...,wn) E H n \ e~ Wl + e~ W2 + ... + e~ WN = 1}. In this case, we 
have Wj(t n ) > for all j. Also there are two indices i,j so that 

(4.4) lim Wi(t n ) = oo and Wj(t n ) remains bounded. 

Now by (4.3) the integral F(w(t n )) is the sum of finitely many integrals of type 
(2.1) over each triangle AviVjVk where the vector (a{ k , af , a^ J ) satisfies the condition in 
corollary 2.3(b). Thus, by the same argument used in the proof of corollary 2.3(b), the 
integral \a{ k — 6j k )dwi + (af — 6f)dwj + (a^ — 6^)dwk is either bounded or tends 
to infinity depending on limsup n ^ 00 max{ r ,s}c{i,j,k}(\ Wr (^ n ^ ~ w s(t n )\) is finite or infinite. 
However, the infinite case must occur due to (4.4). This shows that F(w(t n )) tends to 
infinity. QED 

§5. A Proof of theorem 1.4. 

By the assumption, the solution u(t) = (ui(t), ujsf(t)) of the normalized combi- 
natorial Yamabe flow exists for all time so that there are no singularities forming at 
time equal to infinity. This means that ^(fj's are in some compact interval in R>o 
and also all inner angles 0\ 3 (t) are in some compact interval inside the interval (0, it). By 
corollary 2.2, this implies that the there is a positive constant A so that the eigenvalues 
of coefficient matrix [c rs ]Ar X 7v considered as a bilinear form restricted to the subspace 
{w G R-^ltui + ... + wn =0} is always bounded by —A for all time t G [0, oo), i.e., 

N 

(5.1) CjjWjWj < — A^^iuf, when ^^Wi = 0. 

i,j i i=l 
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To prove the theorem, it suffices to show that the curvature Ki(t) — K av converges to 
exponentiaiiy fast, i.e., there is a positive constants c\, c<i so that 

(5.2) \Ki(t)-K av \ <c ie - C2t . 

Assuming (5.2) holds, then we can solve Ui(t) from the combinatorial Yamabe flow 
and get 

/,\ - / (Ki(s) — K av )ds 

Ui(t) = e J o 

This shows that lim^oo Ui (t) is a positive real number for all indices. Thus the metrics 
u(t)d converges to the constant PL curvature metric. 
To establish (5.1), let us consider the function 

N 

Its derivative can be calculated as 

(5.3) G'(t) =2'£ i c ij (K i - K av )(K J - K av ) 

i,3 

By (5.1), we have G'(t) < -\G(t). Thus G(t) < Ce~ xt . This establishes (5.2). QED 

5.2. Remark. One can now show that the normalized combinatorial Yamabe flow 
for a single Euclidean triangle = — (re/ 3 — dijUi converges to the equilateral triangle 
as follows. By the same argument as in the proof of theorem 1.3, we see that essential 
singularity never occur. If a removable singularity occurs at time t = L, then the three 
inner angles tend to 0, 0, rx. On the other hand the sum of the square of the curvature 
Y2 r (Ki — K av ) 2 achieves its supremum value only when the inner angles are 0, 0, rx. Thus 
by theorem 1.2(a) that the sum of the square of the curvature is decreasing, removable 
singularities never occur. By theorem 1.4, we see the solution converges. This ends the 
proof. 
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Appendix A. A Proof of Theorem 2.1 

We will carry out the computational aspect of the proof of theorem 2.1. 

The following lemma was established in [CL], Lemma A-l. 

Lemma A-l. Suppose AviVjVk is a Euclidean triangle of area A so that the inner angle 
at Vi is 9i and the length of the edge VjVk is Xi. Then 9i is a smooth function of (x^, Xj.Xk) ■ 
The partial derivatives of the function are given by, 
(a)^ = Xi /(2A). 

^H = -gcos(0 fc ). 

In our case, we fix a set of positive numbers di, <ij, dk and choose a conformal factor 
(ui,Uj,Uk) G R->o- The edge lengths of the triangle are Xi = diUjUk- Thus = and 

dx ' i 
dv~ = X i/ U 3- 

Now the partial derivative can be calculated by the chain rule, 

90i ^ d6i dx s 
du r ^ dx s du r 

_ d®j x s 
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This shows that j^-u r = J2 s ^r T^T x s- Now use the l emma A-l above together with 
the fact that Xi = Xj cos 9 k + %k cos 9j, we obtain §^u r = — |^ as stated in theorem 2.1. 

Appendix B. An Application of the Feasible Flow Theorem to Theorem 1.1 

We now verify that the condition (*) in theorem 1.1 is also sufficient for the existence 
of constant curvature PL metric associated to the triangulation. This proof follows the 
same ideas appeared in [Cv]. 

Let us begin with the feasible flow theorem for network flow. Suppose G = (V, E) 
is a directed graph, i.e., a graph so that each edge is oriented. Here V is the set of all 
vertices and E is the set of all oriented edges. For each oriented edge x G E, let int(x) 
and end(x) be the initial vertex and the end vertex of the edge x. (It is possible that 
int(x) = end(x).) For any subset I C V, let int(I) = {i£ E\int(x) G /, end(x) I}, and 
let end(I) = {i£ E\end(x) G I,int(x) I}. 

Assume there is a lower capacity bound a : E — > [— oo, oo] and an upper capacity 
bound b : E — > [— oo, oo] so that a(x) < b(x) for all x G E. 

A feasible flow on the graph G is a function : E — > (— oo, oo) so that Kirchoff's 
current law is satisfied, i.e., for each vertex v 

xEend(v) xEint(v) 

and a(x) < <j>(x) < b(x) for all x E E. 

Feasible Flow Theorem. A feasible flow exists if and only if for every non-empty 
subset U C V so that U ^ V, 

(bi) y, b ( x ^ E 

xEend(U) xEint(U) 

We now apply the theorem to show that condition (*) in theorem 1.1 is a necessary 
condition. To this end, let the set of all vertices in the triangulation be V and the set of 
all triangle in T be F. If v G V and / G F, then f > v means v is a vertex of /. Let z 
be an extra vertex. Define an oriented graph G as follows. The set of all vertices in G, 
denoted by G° is V U F U {z}. The set of all oriented edges G 1 is {(/, v) \f G F, v G V, f > 
v}U{(zJ)\f EF}U{(v,z)\vEV}. 

Now suppose there is a constant PL flat metric associated to the triangulation. Define 
a feasible flow <fi : G 1 — > [0, oo) on the graph G as follows: </>((/, u)) is the inner angle of the 
triangle / at the vertex v, and <j>((z, /)) = 7r and 0((u, z)) = 7r|F|/|V|. Define the lower 
capacity a for G as follows: a((/, u)) = e > which is smaller than any of the inner angles, 
a(x) = (f>(x) for all other edges. Define the upper capacity b for G as follows: b((f, v)) = oo, 
b(x) = a(x) for all other edges. 

To verify the Kirchoff's law for 0, it suffices to check the following three statements. 
At each vertex v G V, it states that the sum of the inner angles at the vertex is 2n — 
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27r(x(M)/\V\) = 7r|F|/|V|, i.e., the curvature is constant. At each vertex / G F, it states 
that the sum of the inner angles of a triangle is n, and at the vertex z, it is the trivial 
statement that n\F\ = |V|(7r|F|/|V|). 

Recall that for a non-empty subset I of V so that I ^ V, we use Fi = {f E F\f > v 
for some v G J}. For this set 7, we consider the proper subset £/ = I U F 7 of vertices 
G°. We claim that the feasibility condition (Bl) for this U is exactly the condition (*) in 
theorem 1.1. 

Indeed, we have end(U) = {(f,v)\v el,f £F r ,f> v}L){(zJ)\f G F 7 }. Note that by 
the choice of F 7 , the set {(f,v)\v G /, / ^ F 7 , / > v} = 0. Thus end(*7) = {(^,/)|/ G F 7 }. 
Also, we have int(U) = {(v, z)\v G 1} U {(/, G F 7 , v I, f > v}. Thus the feasibility 
condition states that 

J]7r>^7r(|F|/|y|)+ e>5>(|F|/|V|). 

/GFj «el {(/,«)l/eFr,t;^/,/>t;} »€/ 

Thus we have |F 7 |/|/| > |F|/|V|. 

It can be shown easily that the condition (*) in theorem 1.1 is exactly the feasibility 
condition in the feasible flow theorem. Thus one can in fact avoid using the result from 
[CL] in the proof of theorem 1.3. Namely, we choose the set of inner angles a{ k in the 
proof by applying the feasible flow theorem on the same graph G and the same upper and 
lower capacities a and b. 
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